Electronically driven superconductor insulator transition in electrostatically doped 

La2Cu04+5 thin films 



J. Garcia-Barriocanal^^^. A. Kobrinskii^. X. Leng^, J. Kinney^. B. Yang^, S. Snyder-*^ and A. M Goldman-'^ 
^School of Physics and Astronomy, University of Minnesota, Minneapolis, Minnesota 55455, USA and 
^GFMC, Dpto. de Fisica Aplicada III, Universidad Complutense de Madrid, 28040 Madrid, Spain 

Abstract 

Using an electronic double layer transistor we have systematically studied the superconductor to insulator transition in 
La2Cu04+i thin films growth by ozone assisted molecular beam epitaxy. We have confirmed the high crystalline quality of the 
cuprate films and have demonstrated the suitability of the electronic double layer technique to continuously vary the charge 
density in a system that is otherwise characterized by the presence of miscibility gaps. The transport and magneto-transport 
results highlight the role of electron-electron interactions in the mechanism of the transition due to the proximity of the Mott 
insulating state. 



I. INTRODUCTION 

Despite extensive efforts, a generally accepted theory 
of the superconductivity of cuprate superconductors is 
still lacking.EI^l The description of the superconductor 
to insulator transition (SIT) has not been completed 
for these materials, leaving important questions concern- 
ing the origin of the High Temperature Superconducting 
(HTS) state, and the phases that are found in its in- 
triguing phase diagram, unresolved.'^ The superconduct- 
ing condensate in cuprates emerges below a critical tem- 
perature (Tc) with the addition of charge carriers to an 
antiferromagnetic Mott Insulator (MI) . It is generally ac- 
cepted that in the vicinity of the Mott state the essential 
2D character of the superconducting Cu02 planes gives 
rise to new properties that cannot be explained by Fermi 
liquid theory.!^ Interestingly, this unconventional behav- 
ior is not limited to the superconducting low temperature 
phase of the cuprates. It also affects the "normal" high 
temperature regime denoted as the pseudogap (PG). Un- 
derstanding the mechanisms that turn Mis into HTSs, 
specifically in the PG region, remains one of the great 
intellectual challenges of condensed matter physics. 

The SIT is one example of a quantum phase transition 
(QPT), which occurs at zero temperature.^ In disordered 
ultrathin metal films it has been interpreted with either a 
bosonic or a femionic description and the consequences of 
each picture are different. On the one hand, the bosonic 
picture involves the presence of Cooper pairs (i.e. the 
charge carriers of the superconducting state) and vortices 
(their resistive counterparts) in both the insulating and 
superconducting phases.^ The Cooper pairs are localized 
and the vortices are delocalized in the insulating regime, 
and the Cooper pairs condense into a charged superfluid 
on the superconducting side of the transition while vor- 
tices become localized. On the other hand, the fermionic 
picture implies the destruction of superconductivity at 
the SIT by means of pair breaking. Accordingly, it can 
be expected that the charge carriers on the insulating side 
of the transition are single electrons or holes, and thus, 
the fermionic picture would highlight the importance of 



electronic correlations possibly due to the proximity of 
the Mott insulator.^ 

Here we report the systematic study of the SIT of a 
four unit cell thick ultrathin film of La2Cu04+5 ((5-LCO) 
using electrostatically induced charge as the tuning pa- 
rameter of the transition. This has been facilitated by 
the combination of the Ozone assisted Molecular Beam 
Epitaxy (OMBE) growth technique together with the 
fabrication of an Electronic Double Layer (EDL) device. 
OMBE provides high crystalline quality and flat ultra- 
thin films of a robust (5-LCO superconductor (Tc ^ 45 
K as measured at the onset of the transition) grown under 
compressive epitaxial strain on top of SrLaA104 (SLAO) 
substrates. The EDL technique employing ionic liquids 
(ILs) has been used to successfully induce levels of dop- 
ing of the order of 10"'^^ cm~^ IMini and it allows for the 
systematic and exhaustive study of the SIT. This ap- 
proach avoids crystal inhomogeneity due to phase sep- 
aration and does not change the random Coulomb po- 
tential at the different levels of doping. In this regard, 
the work presented here would address some of the fun- 
damental issues of HTS. Among other results, in-plane 
low temperature anisotropy of the electronic properties 
of the film has been found. This anisotropy is a minimum 
at the SIT, which occurs at the same hole concentration 
as a maximum in the Hall resistance measured at high 
temperatures, suggesting that the SIT of electrostatically 
doped 5-LCO films is electronically driven. However, it 
should be noted that the system under study differs from 
the bulk due to its low dimensionality, the presence of 
high electric fields at its surface, and the very different 
method of doping. As such it could become a source of 
interesting new physical phenomena. 



II. DELTA DOPED LANTHANUM CUPRATE {S- 
LCO) 

Stoichiometric La2Cu04 (LCO) is an antiferromag- 
netic Mott insulator and is one of the most exten- 
sively studied precursors of HTS (as well as the HTS 



1 



families derived from it) j ^ l ^^ l ^^ l Among those families 
of doped LCO cuprates, the oxygen doped compounds 
(La2Cu04+5 [J-LCO]) stand out for their annealed dis- 
order vs. the quench disorder that arises in cation sub- 
stituted compounds. The interstitial oxygens (i-0) in 6- 
LCO compounds are located in the spacer La202+5 lay- 
ers that are intercalated between the Cu02 layers. These 
added oxygens are mobile down to relatively low temper- 
atures ( > 200 K) and they are able to find an equilibrium 
arrangement in the crystallographic structure.E^KISI j^ter- 
stitial oxygens are located at the (3, 3) crystallographic 
positions and they can form an "ordered sublattice" un- 
der appropriate conditions of synthesis and treatment of 
the crystal. They exhibit staging order in the direction 
perpendicular to the CUO2 planes, and they form in- 
plane stripes along the main diagonal of the tetragonal 
unit cell breaking the symmetry of the Cu-0 octahedra 
along the main in-plane directions of the pseudocubic 
unit cell.^^ The expected pattern of disorder in this fam- 
ily of compounds is therefore weaker than that obtained 
for chemically substituted compounds.'^ (5-LCO crystals 
are believed to be the simplest cuprate HTS since the i-0 
are not randomly distributed, and moreover, they have 
the highest Tc since the i-0 ar e loca ted far away from the 
superconducting Cu02 planes.fl^IIZl Unfortunately, in the 
study of the phase diagram of (5-LCO bulk compounds 
it is impossible to produce a series of doped samples in 
which the concentration of oxygen is systematically var- 
ied since the system exhibits phase separation over wide 
ranges of oxygen doping.-^^ 



III. EPITAXIAL GROWTH AND STRUCTURAL 
CHARACTERIZATION OF 5-LCO THIN FILMS 

A. Ozone assisted molecular beam epiteixy 

Epitaxial (5-LCO thin films have been grown by means 
of the OMBE technique. Pure ozone is distilled from a 
mixture of O2/O3 gas which is produced with a silent 
discharge ozone generator. The ozone gas is supplied to 
the chamber at a constant pressure value of 3 x 10^^ 
Torr during growth. Further details about the MBE sys- 
tem and the sample preparation can be found in previous 
work.BSlThe samples of the present study were grown on 
(001) oriented SrLaA104 (SLAO) substrates at 750 °C 
with a shuttered growth technique. Samples were slowly 
cooled, over 45 minutes, down to 50 -C in an ozone atmo- 
sphere, to prevent de-oxygenation. In order to character- 
ize the growth mechanism, test samples were also grown 
on (001) oriented SrTiOs (STO) substrates. The results 
from the in situ control of the growth and the structural 
characterization experiments were the same in terms of 
crystal quality, for films grown on top of both SLAO and 
STO substrates. 

The growth sequence starts with the deposition of a La 
layer followed by a Cu layer. In situ control of the layer- 
by-layer growth is provided by Reflection High Energy 
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Figure 1: (Color online) The top image is an example of the 
RHEED pattern obtained after the growth of a 4 unit cell 
thick LCO thin film. The bottom graph shows the oscillation 
of the integrated intensity of the specular rod (red line) dur- 
ing the growth of a 12.5 unit cell sample as a function of the 
number of repetitions (odd numbers correspond to the open- 
ing of the Cu shutter and even numbers to the La shutter) . 
The top graph shows the details of the 4th and 5th unit cells 
of the same growth and it represents the single atomic plane 
control of the growth (see text). The (5-LCO crystallographic 
structure is sketched on top of this panel to better visualize 
the growth sequence. 



Electron Diffraction (RHEED) patterns acquired during 
the growth process (top image of Fig. 1). After the depo- 
sition of the La layer a three-dimensional (3D) RHEED 
pattern, characteristic of a rough surface, can be ob- 
served. Once the deposition of the Cu layer is finished, 
the characteristic scattering rods of a two-dimensional 
(2D) growth are revealed indicating the presence of a 
flat surface (see top image of Fig. 1). In the bottom 
graph of Fig. 1 we show the RHEED oscillations ob- 
tained from the integrated intensity over the area of the 
specular rod during the growth of a 12.5 unit cell thick 
sample. Note that the time units of the X axis of this 
graph have been scaled to the number of times that the 
La and Cu shutters were opened. The count starts at 
zero which corresponds to the first La deposition, and 
so, the odd (even) numbers in the X axis of Fig. 1 cor- 
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respond to the closing (opening) of the La shutter and 
the opening (closing) of the Cu shutter. The RHEED 
oscillations reveal that the growth method succeeds in 
preserving a smooth sample surface up to the thickest 
(^300 A) films that were grown. Moreover, one can fol- 
low the growth by counting the number of single atomic 
planes of each element deposited. The middle graph of 
Fig. 1 shows a detail of the same growth corresponding 
to the deposition of the 4th and 5th unit cells of the film. 
The black line in the figure is obtained by integrating the 
intensity of the small spot that is characteristic of the 3D 
growth of the RHEED pattern during the La and Cu de- 
position. It reveals the formation of two LaO planes that 
takes place before the completion of half of a single unit 
cell. The picture shown on top of this graph is a guide 
to following the growth sequence and it is a sketch of the 
crystallographic structure of LCO. 



B. Structural characterization of the LCO thin 
films 

The structural characterization of the films has been 
carried out by means of X-Ray Diffraction and Atomic 
Force Microscopy (AFM). In Fig. 2 we show the X-Ray 
results obtained for a set of characteristic samples. Fig- 
ure 2A shows the X-Ray Refiectivity (XRR) patterns of 
12.5, 12, 4 and 4 unit cell thick LCO films grown on 
STO, SLAO, STO and SLAO substrates respectively, 
from top to bottom. Finite size thickness oscillations 
corresponding to the total thickness of the film can be 
clearly observed over a wide range of values of 29. The 
presence of interference fringes indicates that both the 
substrate-sample and sample-air interfaces are fiat and 
smooth. A quantitative analysis of the XRR results was 
carried out using GenX software. The refinement re- 
sults, the red lines in Fig. 2A are in agreement with 
the nominal thicknesses observed during the growth us- 
ing the RHEED technique. We obtain a roughness of the 
sample-substrate interface of less than 2 A, as well as a 
sample-air interface roughness of about 6 A which are in 
good agreement with the height of half a unit cell. Both 
values of roughness are the same for all of the samples. 

The Wide Angle X-Ray Diffraction data (Fig. 2B) 
show that the sample is grown textured in the (0 1) 
crystallographic direction of the substrate with high crys- 
talline quality. (From top to bottom we show the data 
for the 12.5 and 12 unit cell samples grown on STO and 
SLAO respectively.) No peaks of secondary phases can 
be seen over the whole scan (12 < 2 9 < 65) but peaks 
corresponding to the (0 n) family of crystallographic 
planes of both the substrates and (5-LCO thin films are 
observed. As marked with red arrows in Fig. 2B, the 
position of the (0 4) Bragg peak of the thin film de- 
pends strongly on the substrate used for the growth. As 
expected for growth on a SLAO substrate, an in-plane 
compressively strained thin film with a c-axis lattice pa- 
rameter c = 13.3 A is obtained. 
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Figure 2: (Color online) (a) XRR data of four different 5- 
LCO thin films, from top to bottom: 12.5 unit cells grown 
on STO, 12 unit cells grown on SLAO, 4 unit cells grown 
on STO and 4 unit cells grown on SLAO. Both STO and 
SLAO are (0 1) oriented substrates. The red lines are the 
refinement of the XRR data, see the text for the results of 
the refinement, (b) Wide Angle X Ray Diffraction of the 12 
and 12.5 unit cells samples. In the inset the rocking curve 
of the first Bragg peak of the film (red arrows in figure b) of 
the 12 unit cell sample is shown as a representative example 
of all the samples of this work. The rocking curve has a full 
width at half maximum of 0.03 degrees, corresponding to a 
low crystallographic mosaicity with large domains. 
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Figure 3: (Color online) Reciprocal space maps of the 4 unit 
cell (5-LCO thin film along the (0 1 7) and (1 7) crystallo- 
graphic directions. 



The presence of ozone and the compressive strain re- 
sult in the formation of a (5-LCO thin film where the unit 
cell volume is not preserved and a non-stoichiometric su- 
perconducting phase with excess of oxygen is stabilized 
instead of the insulating antiferromagnetic stoichiometric 
LCO phase. The introduction of oxygen interstitials at 
the La202+5 spacer layers together with the compressive 
strain produces an expansion of the unit cell in the out- 
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Figure 4: (Color online) (Top image) Surface height AFM 
scan of the 4 unit cell thick film (5.3 nm). (Bottom graph) 
Cross-sectional analysis of the diagonal cuts of the top image 



of-plane direction and turns the sample into a high Tc 
superconductor. The samples grown on STO are tensile 
strained and have a c-axis parameter of 13.2 A. With the 
samples grown on STO we checked that the oxidizing at- 
mosphere by itself was able to stabilize (5-LCO regardless 
of the strain state induced in the thin film. In the inset 
of Fig. 2B we show the rocking curve corresponding to 
the (0 4) Bragg peak of the 12 unit cell thick film as a 
representative example. The full width at half maximum 
(FWHM) of the peak is 0.03 degrees which indicates the 
very low mosaicity of our samples and the structural co- 
herence of the thin film in the direction of the growth. 

We use reciprocal space maps (RSM) to get access to 
the in-plane lattice parameters of the sample. In order to 
check both the a and h lattice parameters, we map the 
off-specular (10 7) and (0 1 7) H and K directions of the 
RS respectively. The results for a 4 unit cell thick film 
grown on SLAO are exhibited in Fig. 3, which shows the 
out-plane (L) direction of the RS plotted vs. the in-plane 
one. (Fig. 3A corresponds to the K direction and Fig. 
3B corresponds to the H direction.) The intense peak in 
both maps corresponds to the SLAO substrate refiection 
while the thin film reflection is the wide (narrow) maxi- 
mum along the Y axis (X axis). The maps confirm that 
the sample is epitaxial and are evidence of an in-plane 
coherence length of the thin film limited by the substrate 
in-plane coherence length. The wide thin film peak in the 
L direction is a consequence of the finite size thickness 
(~ 5.5 nm) of the sample. 

Figure 4A shows a 1 /im x 1 /im AFM scan of the 4 
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Figure 5: (Color onhne) (a) Sketch of the EDLT device. The 
top and bottom pictures are the top and side views of the 
device respectively, (b) Resistance versus temperature curves 
of the LCO thin film during the device fabrication and opti- 
mization process. 



unit cell film. The average root mean square roughness, 
arms is 0.2 nm, which is strictly the same as that obtained 
for the 10 /zm x 10 /im scans. The thickness variation 
relative to the total thickness of the sample is 3.75 %. In 
Fig. 4B the height profile scans obtained from the main 
diagonals of the image can be observed. They reveal that 
the largest height variation is ^ 4 A, which corresponds 
to a third of a (5-LCO unit cell. The surface roughness of 
the thin film matches the surface roughness of the SLAO 
substrate after it is annealed under conditions equivalent 
to those of the growth. 



IV. ELECTRONIC DOUBLE LAYER TRANSIS- 
TOR FABRICATION AND OPERATION 

All the set of transport experiments described in this 
work were carried out using the same EDLT. The device 
was fabricated using a 4 unit cell (^5 nm) thick film of 
(5-LCO. To preserve the sample surface during the fab- 
rication of the device, the cuprate film was kept in an 
inert N2 gas atmosphere inside a small portable vacuum 
chamber. When procedures open to the atmosphere were 
required, a continuous fiow of the same gas was used to 
minimize the exposure of the surface to moisture, CO2 
and other atmospheric gases that would potentially de- 
grade the surface quality of the thin film. Silver elec- 
trodes were deposited on the four corners of the thin film 
to allow us characterize its electrical transport properties. 
The electrodes were deposited by a thermal evaporation 
method and during this process the sample was kept at 
liquid N2 temperature to avoid the de-oxygenation of the 
thin film. As sketched Fig. 5a, the sample with the de- 
posited electrodes was mounted in a glass cylinder that 
contains the DEME-TFSI ionic liquid used as the gate 
dielectric. Electrical wires were attached to the Ag elec- 
trodes with indium dots and a Pt coil (gate electrode) 
was installed on the top of the device. 
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In Fig. 5b we show the evolution of the R(T) curves of 
LCO film during the fabrication and optimization of the 
device. Before applying the ionic liquid to the surface 
of the film, and after the above described manipulation, 
the as-grown sample has a R(T) curve characteristic of a 
high quality (5-LCO crystal. The superconducting tran- 
sition temperature (Tc) of the sample is ^-^45 K which is 
the highest Tc reported for thin films of this thicknes^^Hl 
or for the bulk form of this material.l^H Once the ionic 
liquid is added to the cylinder glass container we again 
measured R(T) of the film (red curve). There is a drop 
in the onset temperature of the superconducting transi- 
tion of about 10 K and the transition becomes broader. 
The ionic liquid is added at room temperature and we 
cannot avoid what appears to be a slow room tempera- 
ture chemical degradation of the sample surface, which is 
completely arrested when the sample is cooled down to 
the gating temperatures (240-245 K). We start electro- 
statically charging the (5-LCO film at 240 K by applying 
positive gate voltages (Vg) that induce negative charges 
in the film, and therefore, reduce its carrier concentra- 
tion. The gating process at this temperature and up to 
1.2 V does not produce any significant change in the nor- 
mal state resistance and does not affect Tc . In order to 
build up more charge on the sample surface and enhanced 
the charge transfer to the sample, we increased the ionic 
mobility by increasing the gating temperature to 245 K. 
After recovering the previous Vq, i.e. 1.2 V, the normal 
state resistance of the sample had increased by a factor 
of 1.5, and Tc dropped more than 1 K (the blue line in 
Fig. 5B). Subsequently, all the different values of Vg or 
hole concentration presented in this work were obtained 
after gating the sample at 245 K over a period of 10-15 
min. 

In order to monitor the electrostatic charging of the 
sample, we measured the charge leakage current between 
the Pt coil and the sample during the gating process with 
the configuration illustrated in the bottom panel of Fig. 
6a. A Keithley 6517 electrometer was used as both the 
in-series voltage source and as a picoammeter. An ex- 
ample of the leakage current vs. time measurements is 
shown in Fig. 6a. After removing the constant back- 
ground which appears for the longer times of the mea- 
surement, we integrated the curves as a function of time 
to estimate the amount of charge transferred to the sam- 
ple in each charging step (see the inset of Fig. 6a) . This 
allows us to calculate the total accumulated charge for 
each increment of gate voltage. These are the red sym- 
bols in Fig. 6b. Finally, the concentration of holes at the 
(5-LCO thin film for each Vq (black symbols of Fig. 6b) 
can be derived by assuming that the initial concentration 
of holes per superconducting oxygen plane (holes/Cu) is 
given by the Hall number of the initial Yq measured at 
190 K. The initial density of holes for the 1.2 V gated 
thin film at 245 K is 0.1215 holes/Cu and this is in good 
agreement with the estimation of the number of holes ob- 
tained from the parabolic relationship between Tq and 
the hole concentration for the same conditions .1221 How- 
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Figure 6: (Color online) (a) Leakage current measurement as 
a function of time. The data was taking during the 2.15 V 
gating process. Inset: Integration of the leakage current mea- 
surement after background subtraction, (b) Accumulation of 
charge during the complete experiment (left y axis and red 
dots) and total density of charge carriers per superconduct- 
ing Cu (right y axis and black dots). 

ever, we would like to clarify that this is not an accurate 
estimation of the absolute value of charge density in the 
film, but it allows us to estimate the relative evolution of 
carrier concentration with Yq- 



V. TRANSPORT AND MAGNETO- 

TRANSPORT MEASUREMENTS. RESULTS 
AND DISCUSSION 

The electrical resistance of the sample was measured 
using a four-terminal configuration employing the follow- 
ing devices: a Keithley 7001 switching system, a Keithley 
6221 current source and a Keithley 2182A nanovoltmeter. 
The temperature of the device and the magnetic field con- 
ditions of the measurement were controlled in a Physical 
Properties Measurement System (PPMS). 

We designed an experiment to access three different 
measurement configurations of the in-plane resistance 
based on a four-terminal transport experiment. The first 
two resistance configurations were obtained by apply- 
ing the electrical current in the main tetragonal in-plane 
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crystallographic directions of the film that coincide with 
the edges of the sample. We denote these resistances 
as Rl and R2, as well as their sheet resistance counter- 
parts as Rsl and Rs2 respectively. In order to determine 
the values of sheet resistance along the main crystallo- 
graphic directions of the (5-LCO thin film (Rsl and Rs2) 
from Rl and R2. we applied the method explained in the 
reference. The method is a revision of the Montgomery 
methodf^Sl and it allows us to determine the 2D electri- 
cal sheet resistance tensor of the tetragonal structure of 
the 5-LCO thin film. The third resistance was the Hall 
Resistance denoted by Rxy. We have characterized the 
three resistances in zero field and in a 9 T applied mag- 
netic field. This measurement configuration allows us to 
obtain the ratio Rsl/Rs2. which will be used to evaluate 
the in-plane anisotropy of the electrical properties. 

Figure 7 shows a false color plot of the dependence 
on temperature of Rsl and Rs2 from 3 K to 200 K for 
different values of Vg for both T and 9 T magnetic 
fields. The different values of Vg range from 1.2 V (pur- 
ple or cold) to 3 V (brown or hot) corresponding to a 
concentration of holes per copper oxygen plane of 0.1215 
holes/Cu and (-0.00625) holes/Cu respectively. Figures 
7a and 7c show the data obtained for Rsl in zero field 
and in a 9 T magnetic field, respectively, and Figs. 7b 
and 7d show the results obtained for Rs2 under the same 
field conditions. We display a total of 75 different values 
of Yq, i.e. doping levels, for the data obtained without 
magnetic field (for both Rsl and Rs2) and 33 values of 
Vg for the sets of data taken at 9 T. The legend displays 
the color code for some representative values of Vq. 

As can be seen in Fig. 7, both curves of Rsl{T) and 
Rs2{T) show a characteristic evolution from supercon- 
ducting to insulating behavior as Yq is increased and 
the charge carrier concentration is depleted. For Vg=1-2 
V, the sample is metallic down to the superconducting 
transition temperature Tg~40 K for the T data and 
down to Tc--25 K for the 9 T data. With a further 
increase of Vg the transition temperature slowly moves 
towards lower temperatures as expected for the depletion 
of holes in the superconducting Cu02 planes. Below 40 
K (25 K) the transition is continuous and the correspond- 
ing gate voltage that separates the superconducting and 
insulating regimes is Vg=2.64 V (Vg=2.54 V) for the 
measurements obtained with a T (9 T) magnetic field. 
Under T and 9 T conditions of applied magnetic field, 
Rsl and Rs2 show a similar evolution with Vg- 



A. Finite size scaling 

These data of resistance vs. temperature at different 
gate voltages resemble that of the reported continuous 
SITs as a function of charge carrier concentration in the 
cuprates YBa2Cu307_5^^ and La2_xSrxCu04.^® Con- 
tinuous SITs are quantum phase transitions (QPTs) with 
quantum critical points separating the superconducting 
and insulating ground states. Near the QCP and at 



nonzero temperature the different phases are separated 
by a Quantum Critical Regime (QCR). In the QCR, the 
dependence of a macroscopic observable property, Rs, as 
a function of the tuning parameter, the concentration of 
holes per superconducting Cu02 plane (p) , should follow 
a scaling law of the form: 

Rsip,T)^Rscf{\p-pc\T-^) (1) 

Here Rsc and pc are the critical sheet resistance 
and the critical carrier concentration respectively, T is 
the temperature and vz is the critical exponent prod- 
uct that characterizes the universality class of the phase 
transition.^ 

Following the procedures described in the reference,'^^ 
we first obtained the critical values of carrier concentra- 
tion for the four different configurations of the experi- 
ment. After plotting Rs as a function of p for the lower 
temperature isotherms of the experiment (3-10 K), we 
found a crossover, which is the value of the temperature 
independent sheet resistance at critical doping, pc- (See 
the inset of Fig. 8a as an example.) The value of pc for 
the T measurements is 0.04 holes/Cu and for the 9 T 
measurements is 0.055 holes/Cu. It is independent of the 
in-plane direction of the measurements. We then obtain 
vz by the evaluation of the derivative of Rs with respect 
to p aX Pc- This procedure involves the fitting of all the 
isotherms of the study by a polynomial function that ac- 
curately takes into account the curvature of the data at 
this point. The results of the polynomial fit are the lines 
plotted in the inset of Fig. 8a. We then obtain the criti- 
cal exponent product of the transition as the negative of 
the inverse of the slope of the straight line yielded from 
the log-log plot of (^^)pc ■^s- ^- This is shown in Fig. 
8a along with the linear fits to the data. Interestingly we 
obtained two different values of vz for the T measure- 
ments: vz =1.2 for Rsl and vz =1.4 for Rs2, while the 
transition under a 9 T magnetic field is characterized by 
a single value of vz for both directions of the applied elec- 
tric field: vz = 1.8. As can be observed in Fig. 8a, this 
method is a sensitive tool for determining vz , allowing 
us to distinguish between small variations of the slope. 
Moreover, it allows us to estimate the upper temperature 
limit of the QCR as the temperature at which the anal- 
ysis deviates from the linear regime. Surprisingly, the 
quantum fiuctuations persist up to approximately 30 K. 

From the previous analysis of the data we obtain the 
values of Rsc- The critical value for Rsl / T (Rs2 / 
T) is approximately 21.5 m (15 kn), and for Rsl / 9 T 
(Rs2 / 9 T), RsC decreases to 13.5 kfl (12.5 kQ). The 
measured values of Rsg do not agree with the quantum 
resistance for pairs (6.45 kfl). The result supports the 
possibility of a QPT driven anisotropically in the absence 
of magnetic field and isotropically in its presence. 

An alternative method for determining the critical ex- 
ponent product of the transition is to obtain the best 
collapse of the Rs data inside the QCR onto the scal- 
ing function.'^ In Figs. 8b and 8c we show the results 
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Figure 7: (Color online) Sheet resistance as a function of temperature for different values of Vg and in two different magnetic 
fields. Figs, (a) and (c) show the data taken in the Rsf direction with T and 9 T magnetic fields respectively, and Figs, (b) 
and (d) show the measurements obtain for the Rs2 direction under T and 9 T magnetic fields as well. The inset of Fig. (a) 
shows a magnification of the superconducting part in a log-log scale. 



of this analysis for the T and 9 T transitions respec- 
tively. It highlights the collapse of the data as well as the 
agreement of the values of the critical exponents obtained 
by the two independent methods. This analysis suggests 
that there is a direct SIT between superconducting and 
insulating regimes with the caveat that the analysis is 
cut off at a relatively high temperature leaving open the 
possibility of a more complex scenario at temperatures 
below those actually measured. Indeed, a closer inspec- 
tion of the Rs{T) curves reveals that zero resistance is 
only achieve for the few highest doping levels (see inset 
of Fig. 8a). Otherwise, the low temperature resistance of 
the sample saturates reaching a constant value different 
from the expected O of the superconducting state. This 
result could be a consequence of the occurrence of quan- 
tum fluctuations that would prevent the system from be- 
coming superconducting above K and that would be 
enhanced by the 2D character of the sample .I^Sl 

The SIT in 2D has been the focus of many experimen- 
tal and theoretical works for several years. These include 
studies of amorphous Bi,^® amorphous MoGe,^°'InOx,'^ 



TiN^ STO^ and more recently in LSC0,'2S! YBCd^ 
and FeSe.l231 To our knowledge there has not been any 
previous report of a possible anisotropic QPT that, more- 
over, becomes isotropic under an applied magnetic field. 
This result would imply that there is an intrinsic elec- 
tronic anisotropy in the sample that drives the quantum 
phase transition with a different mechanism depending 
on the in-plane direction. A magnetic field perpendicu- 
lar to the Cu02 planes is able to suppress the anisotropy 
suggesting a link between the superconducting state and 
the anisotropy. To further explore this phenomenon we 
characterized the in-plane anisotropic electronic proper- 
ties of the sample by means of the coefficient Rsl/Rs2. 

B. In plane anisotropy of the sheet resistance 

Figures 9a and 9b show false color contour plots of the 
ratio Rsl/Rs2 as a function of carrier concentration and 
temperature for the data obtained in zero field and in 9 
T. In both maps, cold colors (violet) represent a value of 
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Figure 8: (Color online) (a) Value of the normalized derivative 
of RS as a function of p evaluated at the critical concentra- 
tion of charge density for Rsl (red) and Rs2 (black) in both 
T (solid squares) and 9 T (open circles) magnetic fields. 
Inset: Sheet resistance as a function of charge density for 15 
isotherms ranging from 3 to 10 K (selected temperatures are 
indicated in the legend). In the left panel the sheet resistance 
as a function of the scaling variable is shown for data obtained 
without magnetic field (b) and under a 9 T field, (c). The 
scaling is shown for 6 (12) different hole concentrations rang- 
ing from 0.046-0.034 (0.064-0.034) for data obtained with 
(9) T. The temperatures for both calculations range from 3 
to 10 K. 



Rsl/Rs2=l, with the value of Rsl/Rs2 increasing with 
the increase of the color temperature up to a maximum 
value Rsl/Rs2=5.4 (red.) It is shown that the electri- 
cal resistance of the sample in the absence of magnetic 
field is isotropic in-plane for a wide range of temperatures 
and hole carrier concentrations. For the highest level of 
doping, p'-^0.12 holes/ Cu, anisotropy appears once the 



sample is cooled down to Tc- As the charge density 
decreases, the intensity of the anisotropy increases and 
shifts towards lower temperatures following the line of the 
superconducting dome. At the high carrier concentration 
and low temperature portion of the color map, Rsl /Rs2 
reaches a maximum value of 4 at p=0.06 holes/Cu and 
T=3 K. With further depletion of holes, the intensity of 
the anisotropy falls and it achieves a minimum at the 
critical carrier concentration of holes for the SIT and at 
the lowest temperature of the experiment, i. e., p=0.04 
holes/Cu and T= 3K. From this point downward to the 
lowest carrier concentrations the anisotropy continuously 
increases up to the highest value Rsl/Rs2=5.5, deep in- 
side the insulating regime. In this region of doping the 
anisotropy is present up to 200 K and it decreases as the 
temperature increases. 

After applying a 9 T magnetic field perpendicular to 
the Cu02 planes (Fig. 9b), the anisotropy of the high 
charge density part of the color map is shifted down- 
ward, disappearing from the temperature window of the 
experiment, in a similar manner to the suppression of 
Tc- Interestingly, the rest of the contour plot remains 
strictly the same as the one obtain for T and the only 
differences emerge after cooling the sample below 40 K 
(The coefficient of both Rsl /Rs2 obtained with and 9 
T is shown in Fig. 9c.) 

The presence of in-plane anisotropy in the transport 
properties is a striking result. The anisotropy is devel- 
oped at the superconducting transition which probably 
is enhancing the sensitivity of its detection. Moreover, 
the presence of anisotropy would imply that the macro- 
scopic symmetry of the sample must be broken in order 
to set a preferential channel for conductance along one 
of the Cu-O-Cu directions. In this regard, the struc- 
tural characterization of the sample has not revealed any 
structural or morphological difference along the a and 
b directions of the thin film. Neither AFM nor X-Ray 
experiments suggests this possibility. However, we have 
no tool to characterize the oxygen interstitial order in 
the ultra-thin film and we may be missing an important 
piece of information. 

The observation of a quantum SIT that behaves differ- 
ently depending on the in-plane direction, together with 
the emergence of an in-plane anisotropy in the electri- 
cal properties of the sample, is evidence that the elec- 
tronic ground states of (5-LCO thin films are intrinsically 
anisotropic on both sides of the QCP. Inside the QCR, 
the quantum fiuctuations between both regimes become 
more important. They eventually result in a reduction of 
the electronic anisotropy and the system correspondingly 
becomes more isotropic as it becomes more quantum in 
nature (i.e. Rsl/Rs2 decreases as we approach the QCP.) 

One of the possible ways of interpreting the existence 
of a preferred conductance channel in the system can be 
done in terms of the models of electronic and spin order- 
ings and the experimental observations of stripe forma- 
tion in the 214 family of cuprates.'^^ In this scenario, the 
charge that dopes an antiferromagnetic Mott insulator 
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Figure 9: (Color online) False color contour plots of the coefficient of Rsl/Rs2 obtained with T (a) and 9 T (b). The purple 
color represents a value of Rsl/Rs2~l and the red color represent a value of Rsl/Rs2~5. (c) False color contour plot of the 
coefficient of the T and 9 T Rsl/Rs2 values, (d) Rsl/Rs2 value vs. temperature on different magnetic fields: T (green line), 
-4.5 T (red squares), -9 T (black squares), 4.5 T (red line) and 9 T (black line) for a charge concentration of 0.1215 holes/Cu. 



it is not randomly distributed in the sample. It forms 
ordered clusters of charge and spins promoted by antifer- 
romagnetic correlations and charge density modulations. 
The interlocked patterns can vary from a "classical crys- 
tal (CC)" order of charge on the insulating side, to "quan- 
tum electronic liquid crystal phases (QELC)" or charge 
or spin density waves in the superconducting regime.'^ 
In this regard, neutron scattering experiments on stage 
4 and 6 (5-LCO single crystals have shown that a spin 
density wave state appears at the superconducting tran- 
sition temperature, which is enhanced with a perpendic- 
ular applied magnetic field, contrary to our observations 
on electrical anisotropy in field.l^ In the system under 
study, the 2D character of the ultrathin film could be 
enhancing the essential 2D character of the supercon- 
ducting Cu02 planes, which may result in a ID chain 
clustering of holes. Accordingly, stripe formation could 
produce a lower resistance channel for current flowing 
along the chains of holes than for current flowing across 



them. 

The comparison of data with and without magnetic 
fleld is evidence that the nature of the electronic ground 
states on the insulating and superconducting sides are 
different depending upon the magnetic fleld. On the one 
hand, at the lower levels of doping the anisotropy is more 
stable and the 9 T fleld barely has any effect on its tem- 
perature and charge density dependences. These results 
suggest that the lightly doped insulating state could be 
antiferromagnetic, as might be expected from the bulk 
behavior, and a long-range exchange interaction of spins 
would stabilize robust charge and spin orderings, inde- 
pendent of magnetic fleld. On the other hand, the effect 
of magnetic fleld on the superconducting side reduces the 
intensity of the anisotropy signiflcantly and it only ap- 
pears at the lowest temperatures. This result could be a 
consequence of the reduction of T^. As stated earlier, the 
superconducting transition could be enhancing the sen- 
sitivity of detection and the presence of the anisotropy 
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in other regions of the phase diagram could pass unno- 
ticed. Nevertheless, the effect of magnetic field on the 
anisotropy is evidenced in Fig. 9d. Here we have plot- 
ted the value of Rsl/Rs2 for p=0.1215 holes/Cu mea- 
sured with T (green Hue), 4.5 T (red hne), -4.5 T (red 
squares), 9T (black line) and -9 T (black squares) as a 
function of temperature. The anisotropy is developed 
only in the absence of field and is not correlated with the 
transition temperature, which is shifted systematically to 
lower temperatures with the increasing field. We suggest 
that the anisotropy could arise due to electronic inho- 
mogeneity and the possibility of strong spin-orbit scat- 
tering in a highly 2D system could lead to the magnetic 
response of the anisotropy. However, there is neither a 
specific model for the anisotropy in the conductance of 
the thin film, nor a model for its response to magnetic 
field. 



C. Hall Resistance 

In order to characterize the temperature and doping 
dependence of Rxy of the LCO thin film, we measured 
the Hall voltage (V^) with a 9 T applied magnetic field. 
We did this for both positive and negative currents across 
both of the diagonal contacts of the four-terminal con- 
figuration when the current was applied in the opposite 
diagonal. The same measurements were also carried out 
without magnetic field, to check the zero value (see Fig. 
10a), and with 9 T during the warming of the sample 
from 3 K to 200 K. Accuracy tests under -9 T, -4.5 T and 
4.5 T were measured at selected carrier concentrations as 
a control parameter of the experiment (see Fig. 10b). This 
was also done in resistance vs. magnetic field sweeps. 
The data confirm the linear dependence of Rxy on mag- 
netic field and that the sensitivity of the measurement of 

0.1 n. 

Figure 10b shows the temperature dependence of Rxy 
for different levels of doping ranging from 0.12 holes/Cu 
(violet line) to 0.04 holes/Cu (brown line). Over the 
complete range of carrier concentrations, Rxy increases 
with decreasing temperature up to a maximum value. 
Upon further cooling below the maximum region (T ^ 
15 K) and depending on whether the sample becomes a 
superconductor or an insulator, the value of Rxy keeps 
dropping or eventually increases. It is important to note 
that the change of Rxy{T) to a negative slope at the 
lowest levels of doping and temperatures is correlated 
with non-zero values of Rxy measured without magnetic 
field, and therefore is not meaningful. We have further 
explored the previous result and obtained the maximum 
position for each of the Rxy{T) curves as well as for each 
isotherm. These data are displayed in the inset of Fig. 
10b and show that the Hall resistance has a maximum at 
0.04 holes/Cu at a temperature of ^-^70 K. 

To better visualize this surprising result we show a 
contour color plot of the Rxy measurements in the 
temperature-charge concentration plane (Fig. lOd). 



Note that the purple colors are assigned to Rxy^O ft. As 
a result, Rxy measured with a 9T applied magnetic field 
has a drop-like shape with the maximum values of Rxy 
represented in red and located around the 0.04 holes/Cu- 
70 K point. Moreover, by using this representation we 
can follow the evolution of the maximum of Rxy in the 
T-p plane. The results reveal that the position of the 
maximum as a function of p (the cross symbols in Fig. 
lOd) for the different isotherms is roughly linear with a 
small shift to lower p as the charge density decreases 
(from p=0.04 holes/Cu at the lowest temperatures to 
p=0.035 holes/Cu at 200 K). When the maximum is eval- 
uated as a function of T for the different p values, it is 
a curve shape with a minimum corresponding to p^O.08 
holes/Cu at T-60 K. 

The presence of a maximum in the Hall resistance 
could be revealing an electronic phase transition 
which is a surprising and interesting result. The surprise 
comes from the fact that the maximum value of Rxy is 
found at the critical carrier concentration for the SIT at 
T and the QPT would be a consequence of an electronic 
transition which occurs even at high temperature. At the 
temperatures where the maximum is observed, we have 
found that the longitudinal resistances measured with 
T and 9 T are strictly the same. Hence, the Rxy max- 
imum is evidence of a phenomenon that promotes the 
QPT at T. The interesting implication of this discovery 
is that the SIT in electrostatically doped S-LCO would 
be driven by electronic interactions since it would involve 
a change in the electronic properties of the system when 
it transitions from the insulating to the superconducting 
phases. This scenario is further support by the presence 
of a minimum in the anisotropy of the electrical transport 
properties at the QCP. 



VI. CONCLUSIONS 

The results presented here highlight the role that elec- 
tronic interactions play in the transition from the su- 
perconducting state to the insulating regime. The low 
temperature transport scales down to the lowest temper- 
ature of the experiment, which suggests a possible QPT 
that is anisotropic without magnetic field and isotropic 
in the presence of a 9 T field. The electronic nature of 
this anisotropy is reveal by the measurements of the coef- 
ficient Rsl/Rs2, which interestingly are minimized at the 
critical concentration of charge carriers of the SIT. Fur- 
thermore, at this concentration of holes we have found 
a maximum in the Hall resistance at high temperature 
that could be revealing the electronic nature of the low 
temperature transition. 
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4.5 T (blue line) and 9 T (cyan line). The measurements correspond to a charge concentration of 0.1215 holes/Cu. (b) Rxy as 
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